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TECHNICAL MEMORANDUM I358 

CAICUIATION OF THE SHAPE OF A TWO-DIMENSIONAL 
SUPERSONIC NOZZLE IN CLOSED FORM* 

By Dante Cuneolo 

SUMMARY 


The idea is advanced of making a supersonic nozzle by producing 
one, two, or three successive turns of the whole flow; with the result 
that the wall contoxir can he calctilated exactly hy means of the 
Prandtl-Meyer "Lost Solution. " 


PURPOSE OF THE INVESTIGATION 


The subject matter of this paper is based on the artifice of not 
letting the eacpansion waves emanating from one wall of the nozzle be 
reflected from the opposite side, but of cancelling out their effect by 
compensating compression waves emanating from this opposite wall. In 
this way the difficulties attendant upon the intermeshing of the Mach 
waves are avoided, and it is no longer necessary to integrate the differ- 
ential equation of the hodograph by recasting it as a finite difference 
equation, the solution of which is necessarily approximate. The tech- 
nique illustrated here is nothing more than a quite direct application 
of the Prandtl-Meyer relationship for flow around a sharp corner. On 
the basis of this procedure, the calculation of the coordinates of any 
point on the nozzle contour is independent of that for any other point; 
that is to say, this way of handling the problem eliminates the lack of 
precision usually associated with the tail end of the effiisor in com- 
parison with the better acc\jracy at the beginning sections. 

In addition, what is of the utmost importance, is that it is not 
necessary to determine the characteristics of the flow throughout the 
interior of the nozzle, which leads to a tremendous saving in computa- 
tional labor, with no deterioration in accuracy. 


"Sul Calcolo in Termini Finiti dell’Effusore di una Galleria 
Bidimensionale Supersonica. " L*Aerotecnica, Vol. XXXI, No. 4 , 

15 August 1951 > PP. 225-230. 
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Theoretical Backgrovind 

Let us focus our attention now on the expansion which occiirs when 
a stream turns around a sharp corner (see fig, l) . In this figure 1 
the following notation applies: n is the Mach angle, 'I' is the com- 

plement of the Mach angle, v is the angle through which the flow has 
turned as it progresses through its tvim around the comer, so that 
related to this angle of the flow vector is the angle -d , which is the 
angle between the normal to the wall and the radius vector originating 
at the comer and lying along the points in the flow which have turned 
through the angle v. 

Then, upon making use of well-known theoretical relationships, it 
follows that 


d = V + 


tan Xd = X tan > 1 ^ 


X = 


that is. 


V = d - arctan 


Now let us employ this Prandtl-Meyer relationship in order to 
determine the streamlines of the flow. Let us denote a small lineal 
element of a streamline by the symbol ds (see fig, 2 ) , and then we 
may take dp and pdd to represent the components of this element in 
polar coordinates. 

Thus one may write: 

dp . tan Xd 

= tan V = — 

pdd ^ 


dp tan Xd 



tan Xd 
X 
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and consequently 


log p + — log cos Xd 


= const. 


or 


p(cos Xd) 


X 

X2 


const. 


and upon selection of the value of y as l.UO, one obtains: 


p cos 


6 ± 


P* 


( 2 ) 


One of the walls is thus the straight terminal side of the angle 
of expansion while the other curved wall is defined by means of equa- 
tion (2) as far out as the location where the angle characterizing the 
radius vector has ultimately reached the value dM> at which point the 
desired Mach number will have been attained. From here on out this wall 
is also straight and lies parallel to the terminal side of the comer 
angle. As may be seen from reference to figure 3 the minimum length of 
effusor which it is possible to have is one which terminates right at 
the point where the wall stops cxirving, at which point the very last 
expansion wave has just been included in the process of executing the 
total turn. Of course, the depth of the "throat" or critical cross- 
section is the quantity denoted by p*. 

If one wishes to avoid reliance upon a perfectly sharp corner for 
making the expansion, the two curved walls (see fig. 4 ) given by the 
equations 


6 ^ 

Pi ~ = Pi* 


<6 




6 ^ 

P2 “p ~ P 2 * 




( 3 ) 
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may be employed. Of course, in this case the depth of the throat or 
critical-cross-section is given by p* = - p,,,. 

Naturally, for an effusor conf igiiration such as this, the length 
becomes longer than in the case of the sharp cornered type. The expan- 
sion waves which originate from the number 1 wall are exactly comter- 
acted or 'swallowed up" when they meet the number 2 wall, without pro- 
ducing any reflection. It is just for this reason, even in case one 
wants to give some other arbitrary shape to the number 1 wall, that the 
calculations may be carried out in parametric, but closed-form, provided 
the Mach lines are maintained as straight lines. 

Once it is decided what the shape of the mmiber 1 wall is to be, 
one may find a value of v which corresponds to the value of i3 
defining the direction of the Mach wave with reference to the y-axis. 
This value of V locates a point with coordinates and yi on the 

number 1 wall. Consequently, the point lying on the nimiber 2 wall which 
is marked out by this i9-ray will have the coordinates (see fig. 5) ^ 


X 2 = X]^ + p sin t 3 
y 2 = Yi - P cos ^ 


w 


where 


P = P* 




-6 


Obviously it is necessary to have and y]_ given as functions of 

Let us follow through on the details in the case where the nianber 1 
wall is taken to be circular, with radius = pq. Then the value of 
used as the independent parameter, locates a point with coordinates 
(see fig. 6 ) : 


XI = Pq sin V 






yi = p* + pq(1 - cos v) 


( 5 ) 
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wherein the value of v is given hy: 


V = ■8 - arctan 




and the point lies on this 8 -ray. 

Likewise let us examine more fully the case where the number 1 wall 
is parabolic in shape at the start, with a radiias of curvature equal to 
PO at the apex of the parabola (see fig. 7)* Taking the equation of 
this parabola to be 


y = P* 



2po 


then its slope is given by 


tan V 


dx 


X 

PO 


( 6 ) 


Consequently, since 8 is taken as the independent parameter, one 
obtains 



*1 = PO "tan V 




(7) 


yi = p* 


+ — tan^v 


In those cases where the slope dy/dx has more complicated formula- 
tions (these would be cases devoid of any practical interest for that 
matter, in view of the extreme simplicity of construction exemplified by 
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the contour of the munher 1 wall for the configurations Just examined 
in detail) one may still solve the problem, to whatever degree of 
accuracy is desired, by working with the isocline system: 


V = i8 - arctan 





— = tan V 
dx 


(8) 


y = f(x) J 


which gives the slope of the contour applying at each assigned value of 
the parameter d. Then, after carrying out this construction, one may 
find the corresponding values of the coordinates and y2, by means 

of equations (U) , 


If, however, one insists on having the working section not offset 
from the axis of the throat, it is necessary to apportion the total 
angle through which the stream is turned into three pieces. To be more 
precise, let this total deviation of the flow, corresponding to the 
ultimate Mach number to be atta^ed, be denoted by Vjj. Then the three 
partial turnings of the flow woxold be, respectively, an upward t\nm of 


amount 


V‘, a downward turn of amount 



V *, and once again an 


upward turn of amount — V* = Vy - v”. With this arrangement one 

now has an effusor made up of three expansion regions separated by two 
sections of uniform flow, wherein the Mach numbers attain the values M* 
and M", respectively. 


Now let V be the net angle through which the flow has been turned 
up until the time it has reached a certain location in the effusor; and 
let d be the angle linked to this value of v by means of the Prandtl- 
Meyer relationship, equation (l). In addition, let a be the angle 
swept out by the Mach wave with respect to the y-axis. Then, in the 
first expansion region, we shall have that a = d (see figs. 8 and 9). 
The first Mach wave beginning the second region of expansion will be 
characterized by the angle a = 2i ’ - d * , while within the second region 
it will be true that a = ( 2 ^* - d *) + (d - d*) = d - 2 v The last 
Mach wave ending the second region will be characterized by the 
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angle a = - 2v’, and the first wave beginning the third expansion 

region is then given by a * 2i" - (d" - 2v’) = 2d" - 2v" - d" + 2V’ = d" - Vj^. 
For like reasons we have that, within the third region, the Mach angle 
is given by a = (d" - Vjj) + (d - d") = d - Vjj. In summary, we have that, 
in the three expansion regions, the following hold, respectively 


a » d 

a * d - 2V > 


( 9 ) 


a = d - Vm 


Let us consider what the appropriate relationships are, in the case 
of a t\mnel design such as illustrated in figure 8, wherein the sections 
of tunnel wall effecting the expansions, that is to say, the convex por- 
tions of the wall, consist of circular arcs. For this example the value 
of V* must turn out to be calculated in such a way as to make the pro- 
jection of the broken line ABCDEFGHL upon the y-axis have the 
value 


(ABCDEFGHL) y 



where Ojj and o* are the widths of the effusor at the very end and at 
the throat. Writing this condition in explicit form, we obtain the 
expression 


P* 


Pl(l - cos V’) cos d' - Po cos V + 

/■ A ^ 


cos 


6 ^ 


re 


P2 C08 


ft--) 


6 d" 
cos^ ~ 

% 


cos (d" - 2V‘) - 




2 


( 10 ) 
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where v’ and v” are the net angles throxagh which the flow has been 
turned at the end of the first and at the end of the second expansion 
regions ; i . e . , 


V » = - arctan 


tan X -d ’ 


tan 


Vm 

;*’ = V* + — = - arctan 

2 X 


V- 


(11) 


In fact, one has in particular that: 


(AB)y = Pi 

(BC)y = -Pi COS V* 


(CD)y = - |cd| cos ■d' = -p' cos -d' = - 


cos ■S ' 

6 

cos — - 


(DE)y = -P2 cos v‘ 


(EF)y = P2 cos (EF,y) 


P2 cos - v’j 


(FG)y = p" cos (PG,y) = — - cos (^" - 2 v») 

6 ^ 

cos° — 

fe 


(OH)y « Pj COS - v) 
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In this case where such circular sections of wall are considered 
to effect the e: 35 >ansions^ it is necessary to utilize eqixations (5) 
and (i)-) (when suitably modified for treatment of the flow in the second 
and third expansion regions) for getting the solution, but if this 
circular wall contour is replaced by the curve given by means of equa- 
tions (3)# then this case is handled by setting 


Pl* '■ cos ■ cos (i9' - 2V*) + 


cos 


6 11 
ll6 


COS 


6 11 


-’2* 


cos 


6 r 


cos - 2v') + i-__. cos («" - Vjj) - 


cos 


6 

<6 


Pi* 


cos 


6 1m 


cos - Vj{) = 


OM - 0* 


( 12 ) 


provided it is assumed that 


and 


de| = |bc 1 •= 


Pi* 


cos 


6 11 


JOH 1 = |ef| = Pi" 


Pi* 


cos 


6 il 

\|6 


# 
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In fact, one has in partictilar this time that 


(BD)v ^ cos 

cos6 11 

{6 


angle between ED and y = (ED,y) = d • - 2 v’ 


Po* 

(EG) = cos (BG,y) 

Ad" 

C08° — 

u 


angle between 1 g and y = (EG ,y) = d " - 2V ' 
angle between GH and y = (GH,y) = d” - Vjj 
and 


Pi * 

(HL)v = cos 


COS 


6 % 

<6 


wherein the symbol stands for dj^ - Vj^, 


Nvimerical Applications 

For illustrative pirrposes let us svqppose that one wishes to design 
a tunnel which will reach a Mach number of 


M = 2.5 


by means of an effusor such as depicted in figure 6 , or else with one 
with three kinks as in figiire 8 . 
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Setting down the numerical values available to us, we get the 
following listing: 


= arcos ^ = 66° 25.3* 




M = 




arctan 


= (6 arctan ^ = 105° 32.7* 

{6 46 


Vji = 39° 7.4’ 


P* 




6.5919 


= i /l+j£.V 

p* 6 ) 


2.6368 


Let us compute the value of v* which applies to the type of 
circular arc expansion incorporated into the figure 8 design. In this 
case if we let 


P* * P* - 1 

Pi * 0.5 

P2 = P3 = 1 


Pm - p» 


then, from eqiiation (lO), 


2 


0.8l84. 
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Let the function A(i 3 ') represent the left hand side of equa- 
tion (10) (which is legitimate because V*, v", and d" are all func- 
tions of and we then propose to find the value of d* which makes 


A(d’) = O.818L (13) 


In order to effect the solution of equation (13), the author made 
use of the method of approximation by secants (see fig. lO) . Let us 
assume for example that d^^ ' and ^ 2 ' values of the independent 

variable which when substituted into the expression A(d’) renders a 
result which is too small in the first case and one which is too large 
in the second case in an attempt to find a value of ^ ' which makes 
equation (13) hold true. In addition, let d^ ' be a value of the 

independent variable which is much more accurate, and which can be 
obtained by linear interpolation. Then the two secants Pq P3 and 
P2 P3 will cut the horizontal line A(d') = O.818U in two points, 
which may be labelled Sq and S2. A better approximation to the 
exact solution of equation (I3) is then obtained by selecting d' as 
equal to a value d^ which falls within the interval defined by the 
end-points Sq and S2. After discarding one of the more distant 
points on the curve (in the case illustrated by fig. 10, the point dq* 
would be dropped) the whole process is repeated by working with the 
points P2, P3, and P4 to start with again. 

We shall try out this procedure by first selecting the value of d ' 
as = 60°. Then one finds that: 


d' = 60° 

V = 11 . 90 

V" = 31.5° 
d" = 9^.3° 


A = 0.59 < 0.8184 
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The value of ■ 6 " listed here is obtained hy use of the relationship 


v" = d" -• arctan 

\ 



through use of the method of approximation based on interpolation along 
a tangent. According to this method, if a first approximate solution 
which is too large is determined and denoted by ’ then a more 
accurate approximation to the exact value of d" is given by 


•a 


2 


II 



- v" 

dv” 


= di" 


dd 



cot^ 



where the value of v-j^" is obtained by an analogous formula to the one 
given above for v”; i.e., here one uses the relation 


V 


II 


1 


d 1” - arctan 



From the above derived numerical result it is seen that the guess 
d' = 60° is a too small solution for equation (13). Consequently if 
the value d * = 'JOP is now tried, it tvims out that: 

d' = 70 ° 

V* = 16 . 9 ° 

V” = 36.5° 

d" = 101 . 7 ° 

A = 1 . 47 > 0.8184 
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Thus the trial value = 70 ° is a too large solution for equa- 
tion (13). Linear interpolation then results in the value for i 3 ' 
of 62.6°. 

Consequently, using this value it’ results that; 

i8' =62° 36' 

V’ = 13° 5’ 

V" = 32° 39 * 

= 96° 4 ' 

A = 0.762 < 0.8184 

Now continuing the computation to obtain a better approximation by 
means of the method outlined above (see fig. 10), one gets that 

•a’ = 63° 20» 

V* s 130 26.7' 

V" - 33 ° 0.4 ' 
a" = 96° 36.9 * 

A = 0.8183 

Consequently the solution of equation (13) is a* = 63° 20 *' with 
an error which is less than 1 minute. 

In carrying out the actual construction of the effusor contour it 
is necessary to point out that, if 3 is taken to denote the direction 
of the flow with respect to the x-axis, the following sign relationships 
hold: 


d3 = dv in the first (l) and last (ill) expansion regions 
d 3 = -dV in the central (ll) e35)ansion region 


( 14 ) 



NACA TM 1358 


15 


Now perform the integrations indicated by equations (l4), substi- 
tuting the proper limits, and one obtains 


3 = V 

3 = 2v’ - V 

3 = V - vjj 



(15) 


in the regions I, II, and III, respectively. 

If the coordinates of the points B, E, and H are assumed known, 
then the equations (15) afford the means of determining the coordinates 
of the point (xi,yi) which lies on the convex portion of the wall contour 
and which is related to an arbitrarily selected fixed value of the param- 
eter d . In addition, let the former equations (U) be modified to read 


X2=xi+psina 
72 = yi T P cos a 



(16) 


where the + sign in the second equation holds just for the central II 
region of e^ansion. These expressions, together with the equations (9)> 
allow one to conqjute the coordinates (X2,y2) of the corresponding point 
lying on the concave portions of the wall. 

Working with the illustrative case depicted in figvire 9f let the 
starting data be selected, for example, as 


P* = 1 = d# 

Pi* =0.3 

P2* = Pi* + P* = 1*3 
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Upon substitution of these values into equation (l 2 ) , one will 
find that 


B(d*) = 0 . 8 l 81 ^ 


(17) 


Continuing the computation by an entirely analogous procedure as 
used previously, it is fo\md that 

d* = 63° 5.7» 
v* = 13° 19.6 ' 
v" = 32° 53.3* 
d" = 96° 26.2* 

In order to carry out the actual drawing of the effusor shape (see 
fig* 9) the relationships given as equations (9) are again employed and 
the results applying to the three convex portions of the wall will be 
given by, respectively: 

Xx = XB + PI sin a 
71 = YB - Pi cos a 
xi = XE + PI sin a 
71 = 7 E + PI cos a 
XI = ^ + Pi sin a 
7 l = 7 H “ Pi cos a 

where 

Pi* 

P = 

¥ 



3Q 
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The determination of the concave portions of the wall is carried 
out in an analogoios way, except for the mere replacement of the sub- 
script 1 by the subscript 2 . 


CONCLUSIONS 


The advantages accruing from the methods e^ounded here for design 
of the various kinds of effusor lie essentially in the great degree of 
precision with which it is possible to draw in the wall contour, A 
blemish of the method suggested for design of the effusors possessing 
three successive expansion regions (see figs. 8 and 9 ), which rnxist be 
acknowledged, is that the length is longer than what would be found 
necessary if the effusor had been calculated according to the usual 
hodograph method, although the latter would be less exact (see fig. ll) . 
The lengths of effusor, when non-dimens ionalized by reference to the 
throat depth, have the following magnitiides in the cases exemplified 
by figures 8 and respectively: 


= 12.2 1 

P« 


-L = lk.2 

P* 


> 


(18) 


The same ratio in the case of the effusor represented by figure 11 
has the value 9 * 8 , but it will be cut down to only 4.9 if the effusor 
shown in figure 11 is considered to be merely one half of a complete 
symmetric tunnel. 

This defect in the effusors described here ceases to exist in the 
case of the skewed effusor designs first mentioned in this paper. These 
off-set designs can be fruitfully employed nevertheless where there is 
restricted room for the setup, because it should not entail a very great 
deal of trouble to incorporate a suitable compensating kink in the sub- 
sonic part of the tunnel. 

A promising conqjromise design can be obtained by use of an effusor 
having two expansion regions (with equal and opposite amounts of turning) . 
With this configuration it is clear that the working section will be 
lined up parallel with the axis of the throat section, except that it 
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will be offset laterally (see fig. 12). The length of the effusor and 
its lateral displacement will have the following sizes, respectively 



~ = 1 . 1 ^ 
P* 






(19) 


As is evident, this displacement s/p* is not of formidable size, 
while on the other hand the length of the effusor has been somewhat 
reduced, in comparison with the cases illustrated in figures 8 and 9- 

Despite all that has been said, it is still worthy of note that in 
the case of effusors designed for Mach nmbers which are only slightly 
greater than \anity; that is, for effusors whose lengths do not exceed 
that of the working section, the effusor designed with three expansion 
regions can still be employed to good purpose. Such an eff\isor designed 
to produce a Mach ninnber of M = 1.2 is illustrated in figure I3. 


Translated by R. H. Cramer 

Cornell Aeronautical Laboratory, Inc., 

Buffalo, New York 
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